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We investigate the masses of the light quarks with lattice 
QCD. We show that most of the large dependence on the 
lattice spacing a observed in previous determinations using 
Wilson fermions is removed with the use of an O(a) corrected 
action. In the quenched approximation, we obtain for the 
strange quark MS mass m s (2 GeV) = 95(16) MeV, and for 
the average of the up and down quark masses m;(2 GeV) = 
3.6(6) MeV. Short distance arguments and existing staggered 
fermion calculations make it likely that the answers includ- 
ing the effects of quark loops lie 20% to 40% below this: 
m s (2 GeV) in the range 54-92 MeV, and m;(2 GeV) in the 
range 2.1-3.5 MeV. We argue that almost all lattice deter- 
minations of the light quark masses are consistent with these 
values. These low values are outside the range usually given 
by conventional phenomenology. 

PACS numbers: 14.65. Bt, 12.15.Ff, 12.38. Gc 



Among the most important applications of lattice 
gauge theory to particle physics are the calculations re- 
quired to determine the fundamental parameters of the 
quark sector of the standard model. One of the most 
important of these is the overall scale of the light quark 
masses. It is one of the least well known of the funda- 
mental parameters of the standard model. (Estimates for 
the strange quark mass range from 100 to 300 MeV for 
the MS masses renormalized at a "high" energy scale, 
1 GeV, and for the average light quark mass from 3.5 
to 11.5 MeV §.) It is also one for which lattice meth- 
ods are almost uniquely reliable, unlike quark mass ra- 
tios or the strong coupling constant a g , for which other 
powerful methods exist. Values for quark masses have 
been obtained since almost the beginning of lattice phe- 
nomenology However, improved understanding of 
perturbation theory and finite lattice spacing errors has 
been required to make sense of the various lattice deter- 
minations, which initially ranged over a factor of three. 

Lattice determinations of standard model parameters 
consist of two pieces. Calculations of experimentally 
measurable quantities such as hadron masses are used 
to fix the bare coupling constant and quark masses in 
the lattice Lagrangian. Short distance calculations are 
used to relate the bare parameters in the lattice theory 
to renormalized, running coupling constants and masses, 
such as those of the A/5 scheme. 



Quark masses are most easily obtained in lattice cal- 
culations by matching pseudoscalar meson masses with 
experiment. These are among the easiest lattice calcu- 
lations, having small statistical and finite volume errors. 
Experimental uncertainties are also negligible. Uncer- 
tainties are dominated by truncation of perturbation the- 
ory and discretization errors, and by errors arising from 
the omission of light quark loops (the "quenched" ap- 
proximation) . 

The short distance calculations relating the parameters 
in various regulators may be performed by demanding 
that short distance quantities such as the heavy-quark 
potential or current correlation functions be the same in 
both regulators. It is desirable to do the lattice part 
of such calculations nonperturbatively as much as pos- 
sible, to test for the presence of nonperturbative short 
distance effects and possible poor convergence of pertur- 
bation theory. Such nonperturbative short distance anal- 
ysis for quark masses is currently less advanced than the 
analogous investigations for the strong coupling constant. 

Perturbative relations between the lattice bare mass, 
mo, and the MS mass, m, may be obtained by demand- 
ing that on-shell Green's functions calculated with both 
regulators be equal. Analogous perturbative expressions 
for the renormalization of a s were initially rendered al- 
most useless by sick behavior in the lattice perturbation 
series. In Ref. [|J it was shown that such behavior could 
be understood and mostly eliminated by a mean field the- 
ory resummation of large "tadpole" graphs. The effects 
of such large tadpoles are much less important for quark 
mass renormalizations than for a s B. 

To reduce the effects of such graphs further, the ex- 
pression giving m from mo may be rewritten in terms of 
a mean field improved mass m, 



m(fi) = rh 1 + a s 7o (lnC m — In {a/i)j , 



(1) 



where 70 = 2/tt is the leading quark mass anomalous di- 
mension, and In C m is the result of a one loop calculation. 
Here we use m = mo/ {J (Up) for the mean-field-improved 
bare mass m. The nonperturbative value of the plaquette 
expectation value (Up) is used in the expression for m to 
incorporate an estimate for higher order tadpole graphs. 
The one loop term lnC m is then adjusted to remove the 
one loop part of this expression {/ (Up) = (1 — (7r/3)a s ). 
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FIG. 1. Previous lattice results for the MS masses of the 
light quarks, renormalized at 1 GeV, with the lattice spacing 
set by the rho mass. Lattice spacing dependence is large for 
quenched Wilson fermions (diamonds) and small for quenched 
staggered fermions (filled squares). Results from two-flavor 
staggered fermion QCD (open squares) lie below those from 
quenched approximation staggered fermions by a reasonable 
amount. Data from Ukawa pi. 

In Fig. |l| we show a compilation of previous results 
given by Ukawa (|]. Quenched results obtained with 
staggered fermions are almost cut-off independent for lat- 
tice spacings less than 1 GeV -1 . However, for staggered 
fermions the constant in Eq. (|l|) is C m = 132.9 ||. This 
leads to correction factors of 50-100%, most of which 
is unexplained by mean field theory, casting doubt on 
the reliability of the perturbative relation between the 
staggered-fermion quark mass and the M S quark mass. 

For Wilson fermions, we have C m — 1.67 ||, and thus 
a well-behaved perturbation series. However, the numer- 
ical results for the Wilson action show large cut-off de- 
pendence. They lie far above the results for staggered 
fermions but show a downward trend as the lattice spac- 
ing is reduced. The Wilson fermion action contains an 
error of 0(a), which is absent in the staggered fermion 
action. If the results are extrapolated in a, one obtains 
a result much closer to the results of staggered fermions. 
(See, for example, Ref. ||].) However, remaining sources 
of cut-off dependence are an unknown combination of 
0(a 2 ), 0(a s a), 0(a 2 ), etc. They cannot be estimated or 
removed by simple extrapolation, since we do not have a 
quantitative theory of their functional form. One there- 
fore needs to investigate the removal of the dominant 
0(a) error from the Wilson action. 

A convenient action for doing this has been pro- 
posed by Shcikholeslami and Wohlert ||. Their im- 
proved action incorporates an extra dimension five term 
ipPtivFuvip, the so-called "clover" term. The one-loop 
correction to the coefficient of the clover term is large [[l0| , 
as suggested by mean field theory j| . It is a three-tadpole 
correction and can be approximated by c ~ (t/p)~ 3//4 , 
where the tree level coefficient is normalized to be one. 
For the improved action, C m = 4.72 Thus, Eq. (|l]) 



is still well-behaved. 

We use this action to determine the overall scale of the 
light quark masses. (Or equivalently, the coefficient of mi 
in the expression M 2 = Cmi + . . .. We do not see devia- 
tions from the leading order of this equation, see below.) 
Our lattice spacings range from (the coarse) 1.26 GeV -1 
(at which uncertainties due to perturbation theory are 
starting to approach 50%), down to 0.39 GeV -1 (where 
perturbation theory appears well behaved). We have per- 
formed the calculation at the largest lattice spacing to 
investigate its behavior where it is beginning to break 
down, but we omit it from our final results. The lattice 
spacings have been obtained from the IP- IS splitting of 
the charmonium system, Mh c — (3Mj/^, + M^J/4, for 
which the uncertainties of lattice calculations are partic- 
ularly small and easy to understand. This means that nu- 
merical uncertainties in our results for the quark masses 
arise from a combination of uncertainties in the charmo- 
nium and pion calculations. 

We use improved lattice perturbation theory to convert 
to the MS mass at renormalization scale /i = 2 GeV and 
charmonium splittings to determine the lattice spacing, 
whereas previous determinations typically used bare per- 
turbation theory at scale /j, = 1 GeV and the rho meson 
mass to determine the lattice spacing. Although renor- 
malization at 1 GeV is conventional in nonlattice results, 
renormalizing down to such a low scale introduces addi- 
tional perturbative uncertainty into the results which is 
not present in the underlying lattice results. 

Discussions of our charmonium calculations have ap- 
peared in Ref. [Q. Some technical details and results 
of our pion calculations are given in Table Q, For our 
most significant data point, the improved clover action 
at (3 = 6.1, we have used 100 configurations separated by 
4000 heat bath gauge sweeps. Pion correlation functions 
were calculated using 2x2 correlated fits (fitting two 
states using two operators for the pions). Contributions 
from excited states were checked further on the smaller 
lattices by comparing with lxl and 3x3 fits. Statisti- 
cal errors were calculated using 1000 bootstrap samples. 
Longer descriptions of our analyses for the charmonium 
system and for the light quark masses are in prepara- 
tion If]. 

In Fig. | and in Table we show our results for the light 
quark masses in the quenched approximation. The errors 
shown are statistical only. The diamonds are our results 
for unimproved Wilson fermions. They are consistent 
with the existing determinations (diamonds in Fig. ^). 
The triangles are our results for the mean-field-improved 
clover action. Most of the cut-off dependence has been 
removed. 

Remaining sources of such cut-off dependence could in- 
clude large a 2 corrections to the mass relation, Eq. 
further corrections to the clover coefficient in the pion nu- 
merical calculations, and 0(a 2 ) corrections to the char- 
monium IP-IS splitting. O(a) corrections are expected 
to be negligible for this splitting, but 0(a 2 p 2 ) corrections 
could be larger since quark momenta are larger in char- 
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monium than in pions. We estimate 0(a 2 p 2 ) corrections 
in the charmonium splitting to range from 4% to 20% on 
our three finest lattice spacings. The perturbative one- 
loop result for the coefficient of the 0(a) clover correc- 
tion agrees with the mean field estimate fll(J. However, 
a nonperturbative determination appears indeed to favor 
a further significant correction Jlij . Purely perturbative 
errors in the relation between the lattice and M S masses 
should be of order a 2 ~ 5% at our finest lattice spacing. 
Other smaller uncertainties include finite volume effects, 
which are expected to be a couple of per cent or less, and 
statistical errors, which are 4% and arise mostly from the 
lattice spacing derived from the charmonium system. 
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FIG. 2. Our results for the masses of the light quarks. 
Most of the lattice spacing dependence of unimproved Wil- 
son fermions (diamonds) is removed by the use of an 0(a) 
corrected action (triangles) with a tadpole improved coeffi- 
cient. The lattice spacing is set by the charmonium IP-IS 
splitting. Errors are statistical only. 

We have examined the pseudoscalar meson mass 
squared as a function of the quark mass, which should 
be linear plus small corrections in the small quark mass 
limit. Our numerical data are for quark masses in the 
range 0.4to s to m s . In this range, we find no statistically 
significant evidence for quadratic terms in M 2 vs. m;, 
much less the very large quadratic terms that have been 
postulated to make m u = consistent with experiment. 
Therefore, our results for the ratio of the strange to light 
quark masses agree with lowest order chiral perturbation 
theory: (m s + m;)/(2m/) » M^ /M^ sa 13.6. The val- 
ues for to/ are obtained by linear extrapolation from the 
lowest masses at which we have performed simulations 
down to the physical mass value. 

At present, the uncertainty associated with the remain- 
ing cut-off dependence is the least reliably understood 
uncertainty in the quenched approximation. Pending 
further understanding of this error, we take our result 
at the smallest lattice spacing as the top of our lattice 
spacing error bar. We take a linearly extrapolated re- 
sult (the lower of two plausible extrapolation methods) 
through the three finest lattice spacings as the bottom. 



This gives a range of 0.8 MeV for the cut-off dependence 
uncertainty, and we take the center of this range as our 
continuum limit, quenched approximation result: 



m 4 (2 GeV) = 3.6(6) MeV, 
m s (2 GeV) = 95(16) MeV. 



(2) 
(3) 



The perturbative and cut-off dependence uncertainties 
were added linearly in the total error, since they are re- 
lated. All other uncertainties were added in quadrature. 

Another determination of the strange quark mass with 
an O(a) improved action has been reported Jl5| . This 
determination used a tree-level, rather than a mean-field 
improved, estimate for the clover coefficient. They ob- 
tained m s (2 GeV) = 128(18) MeV. They did not attempt 
to correct for the effects of the remaining lattice spacing 
dependence or the effects of the quenched approximation. 
Most of the discrepancy with our results arises from fact 
that we have used much larger clover coefficients, and 
make an allowance for the fact that we continue to find 
significant cut-off dependence even so. 

In the quenched approximation, QCD couplings run 
slightly incorrectly. The strong coupling constant, for 
example, runs too fast without the effects of light quark 
loops |12] ]. To leading logarithmic accuracy, a s (ir/a) is 
too small by a factor of about (3^ / (3^ , where (3^ and 

Pq are the leading quenched and unquenched p func- 
tions, respectively. This means that the running of the 
quark mass in the perturbative momentum region around 
it /a is too slow, by about the same factor. In Ref. [^6| , 
the ratio of quenched and unquenched quark masses aris- 
ing from the perturbative region was estimated, to lead- 
ing logarithmic accuracy, to be 



m(7r/q)| qu . 
ra(7r/a)| 



unqu. 



a.(7r/a)*<W-W> (4) 
1.15 to 1.20, (5) 



for a s {ir/a) w 1/6 to 1/8. There is, of course, an ad- 
ditional contribution from the nonperturbative region, 
which is unknown. However, a correction due to light 
quark loops of tens of per cent in the downward direc- 
tion from the perturbative region at least would not be 
unexpected. 

Some quenched and unquenched staggered results sum- 
marized in Ref. B are shown in Fig. (Unquenched Wil- 
son fermion calculations appear to be much more difficult 
to perform and harder to interpret.) The unquenched re- 
sults indeed lie below the quenched results by roughly the 
expected amount, and we take them seriously enough to 
use them to estimate the effects of quenching. We argued 
above that quenched staggered quark mass determina- 
tions look good in most ways, but are unreliable because 
of the poor convergence of perturbation theory. How- 
ever, the large corrections cancel out in the ratio of the 
quenched and unquenched determinations, making this 
a useful quantity to examine. To minimize effects due 
to differences in analysis methods, we estimate the ratio 
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from the results of a single group, at similar volumes and 
lattice spacings (about 0.4 GeV- 1 ) 0H, 

and obtain 



to; (1.0 GeV) nf = _ 2.61(9) 
m«(1.0 GeV) nf=2 ~ 2.16(10) 
= 1.21(7). 



(6) 
(7) 



Since there are, in fact, three flavors of light quarks in the 
world and not two, we will use this ratio as a lower bound 
on the actual ratio and use the square (corresponding to 
four light quarks) as an upper bound. 

In summary, after making some plausible cuts, existing 
determinations are reasonably consistent, or have plausi- 
ble explanations for discrepancies. We omit results with 
very small physical volumes (smaller than 1.5 fm) and 
very large lattice spacings (larger than 0.2 fm, or 1.0 
GeV -1 ). We also do not attempt to interpret the re- 
sults with unquenched Wilson fermions, which are in a 
more primitive state than those with staggered fermions. 
Of the remaining determinations, we have shown that 
the cut-off dependence and large size of determinations 
with quenched Wilson fermions arise mostly from the 
well-known 0(a) error. The remaining discrepancy be- 
tween the quenched clover-improved fermion results and 
the quenched staggered fermion results is plausibly at- 
tributed to the apparent poor convergence of staggered 
fermion perturbation theory and the remaining cut-off 
dependence in the improved fermion results. The small 
difference between quenched and unquenched staggered 
fermion results is roughly what is expected. Putting all 
this together, we arrive at the following estimates for the 
light quark masses including effects of light quark loops, 
which we believe are consistent with all known facts: 

• to s (2 GeV) in the range 54-92 MeV, 

• to; (2 GeV) in the range 2.1-3.5 MeV, 

for the MS masses renormalized at 2 GeV. These 
estimates arise from combining our quenched result, 
Eq. (§), with the correction ratio obtained from stag- 
gered fermions, Eq. (Q). Renormalizing down to the 
scale 1 GeV, where conventional mass estimates are often 
quoted, the estimates are raised by 10%, to m s (l GeV) 
in the range 59-101 MeV, and to;(1 GeV) in the range 
2.3-3.9 MeV. 
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TABLE I. Our results for m;(2 GeV), the average of the u 
and d quark masses, renormalized at 2 GeV. c is the coefficient 
of the O(a) improvement operator. 



p 


5.5 


5.7 


5.9 


6.1 


a (GeW 1 ) 


1.26 


0.86 


0.57 


0.39 


volume 


8 3 x 16 


12 3 x 24 


16 3 x 32 


24 3 x 48 


m (c = 0) 


6.31(26) 


5.93(17) 


4.88(18) 


4.62(22) 


m (improved) 


4.75 (19) 


4.41(12) 


3.90(13) 


3.84(18) 


c 


1.69 


1.57 


1.50 


1.40 



4 



